A&A manuscript no. 

(will be inserted by hand later) 



Your thesaurus codes are: 

08.16.4, 08.13.2, 13.09.6, 09.16.1, 09.04.1 



ASTRONOMY 
AND 
ASTROPHYSICS 
1.2.2008 



OO 

m 



> ■ 

o^ : 

OO 
(N ■ 
(N ■ 

o : 

OO , 

ON • 

i!^ ■ 
Oh 



Crystalline silicates in Planetary Nebulae with [WC] 
central stars0 

L.B.F.M. Waters^ ^ D.A. Beintema^, A. A. Zijlstra^, A. de Koter\ F.J. Molster\ J. Bouwrnan^ T. de Jong^ \ 
S.R. Pottasch*, and Th. de Graauw^ 

^ Astronomical Institute Anton Pannekoek', University of Amsterdam, Kruislaan 403, NL-1098 SJ Amsterdam, the Netherlands 
^ SRON Laboratory for Space Research Groningen, P.O. Box 800, NL-9700 AV Groningen, The Netherlands 
^ Department of Physics, UMIST, P.O. Box 88, Manchester M60 IQD, U.K. 

Kapteyn Institute, University of Groningen, P.O. Box 800, NL-9700 AV Groningen, The Netherlands 

received date; accepted date 



Abstract. 

We present ISO-SWS spectroscopy of the cool dusty en- 
velopes surrounding two Planetary Nebulae with [WC] cen- 
tral stars, BD+30 3639 and He 2-113. The X < 15 fj.m region 
is dominated by a rising continuum with prominent emis- 
sion from C-rich dust (PAHs), while the long wavelength part 
shows narrow solid state features from crystalline silicates. 
This demonstrates that the chemical composition of both stars 
changed very recently (less than 1000 years ago). The most 
likely explanation is a thermal pulse at the very end of the 
AGB or shortly after the AGB. The H-rich nature of the C-rich 
dust suggests that the change to C-rich chemistry did not re- 
move all H. The present-day H-poor [WC] nature of the central 
star may be due to extensive mass loss and mixing following 
the late thermal pulse. 
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1. Introduction 

Planetary Nebulae (PNe) are the ionized and photodissociated 
remnants of extensive mass loss which the central star expe- 
rienced when it was a cool Asymptotic Giant Branch (AGB) 
star. PNe therefore provide insight in the mass loss on the AGB 
and beyond. An interesting class of PNe are those with a Wolf- 
Rayet ([WC]) central star. These central stars have strongly 
enhanced C and He, but little or no H in their atmosphere, 
and are the low mass counterparts to the population I Wolf- 
Rayet stars. About 50 [WC] stars are known (Corny & Stasin- 
ska 1995). The properties of these nebulae do not differ sig- 
nificantly from nebulae with 'normal' central stars (Pottasch 
1996). 
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The formation of H-poor central stars is somewhat of a 
puzzle, but is probably related to a thermal pulse either at 
the very end of the AGB or young PN phase (in this paper 
referred to as PAGB pulse; Zijlstra et al. 1991), or when the 
star is already on the cooling track (very late thermal pulse; 
Iben 1984). The H-rich layers may be removed due to efficient 
mixing and subsequent nuclear burning, or by extensive mass 
loss, exposing processed layers to the surface. Examples of ob- 
jects that may have recently experienced a late thermal pulse 
are FG SCe and Sakurai's object (e.g. Paczyiiski 1970; Nakano 
et al 1996; Duerbeck et al. 1996; Blocker & Schonberner 1997). 
It is not clear that these stars will develop [WC] type spectra. 

IR spectroscopy of the dust ejected by [WC] central stars 
shows the well-known family of IR emission bands, usually 
identified as due to Polycyclic Aromatic Hydrocarbons (PAHs) 
(e.g. Cohen et al. 1986) . This confirms the C-rich nature of the 
most recent mass loss episode. In this Letter, we present new 
infrared spectra taken with the Short Wavelength Spectrome- 
ter (SWS) on board of the Infrared Space Observatory (ISO) 
of two well-studied PNe with [WC] central stars, BD-l-30 3639 
and He 2-113. Both nebulae have associated molecular gas 
(Taylor et al. 1990; Bachiller et al. 1991; Gussie & Taylor 1995). 
In Sect. 2, we discuss the observations and data reduction. 
Sect. 3 gives an inventory of the solid state features (both C- 
rich and O-rich). In Sect. 4 we present a first attempt to model 
the dust spectrum of BD-l-30 3639, and we discuss these results 
in the context of post-AGB evolution and formation scenarios 
for [WC] central stars of PNe. 



2. The Observations 

Full scans (AOTOl speed 1 and 2, A 2.4-45 ^m) of BD-l-30 3639 
([WC9]) and He 2-113 ([WCll]) were obtained using the SWS 
(de Graauw et al. 1996) on board of ISO (Kessler et all. 1996), 
as part of a guaranteed time programme on spectroscopy of 
PNe (p.i. D.A. Beintema). The spectra were taken on Nov. 
6, 1996 (BD-h30 3639) and February 4, 1996 (He 2-113), and 
were reduced using version 5.3 of the SWS off-line analysis 
pipeline. All detector signals were inspected for spurious jumps 
and glitches, and in two cases removed manually. The spectra 
of the 12 individual detectors in each SWS AOT band were 
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Fig. 1. 2.4-45 fim ISO-SWS spectra of He 2-113 and BD+30 3639, 
showing the PAH emission features at A < 13 /im and the crystalline 
silicates at A > 20 /im. Tickmarks indicate olivine bands (short thin 
lines; Dorschncr, private communication) and pyroxene (long thick 
lines; Koike ct al. 1993). Emission from [Ne II] at 12.8 /um is seen in 
both nebulae, and from [Ar II] at 6.99 ;um, [S III] at 18.71 ^im, and 
[Si II] at 34.8 /im in BD-l-30 3639. 

sigma-clippod, averaged and rc-binned to a uniform resolution 
of 250. Finally, the different AOT bands wore joined together 
to form a continuous spectrum. For He 2-113 we did not need to 
adjust individual AOT bands, but for BD-l-30 3639 we reduced 
the fluxes measured in bands 3E and 4 (wavelengths of 27.5 
microns and longer) by 17 per cent. This is within the absolute 
flux calibration uncertainties (Schaeidt et al. 1997). Jumps in 
flux can occur at 27.5 and 29.0 microns, where the effective 
aperture changes from 14*27 to 20*27 and 20*33 arcsec). The 
jump seen in BD-l-30 3639 could be caused by the finite extent 
of the nebula (10 arc sec or less in the optical, Harrington et al. 
1997); Cox ct al. (1997, irr preparation) firrd that the extent of 
the neutral envelope is twice that of the ionized nebula. Flux 
jumps can also be due to a small pointing error. The good 
agreement with the IRAS-LRS spectrum (not shown) excludes 
a serious pointing error. The resulting spectra are plotted in 
Fig. 1, and in Fig. 2 we plot continuum subtracted spectra for 
both objects. The location of the continuum was estimated by 
eye. This procedure may underestimate the contribution frorrr 
weak, broad features that are not easily distinguishable from 
the smooth underlying continuum. However narrow (width less 
than a few /^m) features are well preserved. 

3. Solid state features 

Here we discuss the solid state features in our SWS spectra. A 
full description of the atomic and ionic emission lines will be 



given elsewhere. 

The short wavelength part of the spectra of both stars is 
dominated by prominent emission from the well-known family 
of IR emission bands, usually attributed to Polycyclic Aro- 
matic Hydrocarbons (PAHs) (e.g. Leger & Puget 1984). PAHs 
in BD +30 3639 were previously reported by e.g. Witteborn 
et al. (1989), AUamandola et al. (1989), and Cohen et al. 
(1986, 1989). The overall shape of the PAH bands is similar in 
both objects (Fig. 2). The plateau at 11-14 fim is very promi- 
nent in BD-l-30 3639, but does not extend as far to the red 
as in IRAS21282-f5050 (Beintema et al.l996; Molster et al. 
1996). This plateau is attributed to a blend of rrrany bands 
due to deformatiorr rrrodes of large PAH molecules. As in 
IRAS21282+5050, we find prominerrt 12.0 and 12.7 ^J.In fea- 
tures. These latter bands may be attributed to C-H bending 
modes in PAH molecules with 2 or 3 adjacent H atoms at- 
tached (AUamandola et al. 1986). This evidence of H in the 
C-rich dust shell is consistent with a recent determination of 
the photospheric H aburrdance in He 2-113 (7 per cerrt, Leuerr- 
hagen & Hamann 1997), although a quantitative comparison 
of the present-day photospheric H abundance and the H con- 
tent of the C-rich dust shell is difficult to make. The PAHs in 
BD-l-30 3639 suggest that also in that star H was present in 
the C-rich phase. 

In both objects, the spcctrurrr at A > 20 ^im shows weak, 
and narrow solid state features, whose strerrgth is of the order 
of 10 per cent of the local continuum (Fig. 2). These solid state 
features are observed in the spectra of oxygen-rich (proto-) 
planetary nebulae and red AGE stars (Waters et al. 1996; Just- 
tanorrt et al. 1996). Following Waters et al. (1996) we iderrtify 
the features at 19.8, 23.5, 27.5 and 33.8 yum with crystalline 
olivine (Koike et al. 1993; Dorschner, private communication). 
The peak at 40.5 /xm is due to crystalline pyroxenes (Koike et 
al. 1993; Dorschrrer, private communication; Jager et al. 1994). 
The pyroxerres also contribute to the peaks at 23 pm and 33 
fim. The barrd strength ratio of the 23.5 and 33.8 /xm bands 
of BD-l-30 3639 and He 2-113 (0.5 and 1 respectively) suggests 
that the crystalline dust in He 2-113 is warmer than that in 
BD+30 3639 (assuming similar dust composition). This is con- 
sisterrt with the slope of the continuum. We corrclude that both 
[WC] stars experierrced strorrg mass loss as an 0-ricli star at 
the end of the AGB, and that the change to C-rich chemistry 
immediately followed. The 0-rich nature of the outer layers 
of the BD-l-30 3639 nebula can explain the high albedo of the 
scattering particles in the halo surrounding the ionized inner 
nebula (Harrington et al. 1997). 

4. Discussion 

4.1. A dust model 

The presented model is calculated using the dust radiation 
transfer code modust (de Koter et al. in preparation). We 
have modeled the SWS spectrum of BD+30 3639 using a two 
comporrent dust rrrodel, of which the inrrcr and outer shell con- 
nect. The inner part contains amorphous carbon grains and we 
have constrained the radial dimensions of this shell by requir- 
ing that it coincides with the region from which PAH emission 
is observed (Bernard et al. 1994), i.e. from about 1.9" to 2.7". 
Adopting a distance of 2.6 kpc (Hajian & Terzian 1994), this 
yields an inner and outer radius of the C-rich shell of 7.3 10^^ 
and 10.6 10^® cm respectively. The luminosity of BD+30 3639 
was set to 22,000 Lq and its effective temperature to 34.2 kK. 
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Fig. 2. Continuum subtracted 2.4-45 /im spectra of He 2-113 (top) 
and BD+30 3639. Notice the prominent UIR emission bands at 3.3, 
6.2, 7.7, 8.6, and 11.3 /im, often attributed to Polycyclic Aromatic 
Hydrocarbons. At longer wavelengths emission from crystalline sili- 
cates is found (olivines, pyroxenes). 



Both these values are derived from the analysis of Siebenmor- 
gen et al (1994) after scaling to our (larger) adopted distance. 
Note that Tea is rather uncertain as estimates of this value 
range from 30 kK (Siebenmorgen et al. 1994) to 47 kK (Leuen- 
hagen et al. 1996) . The stellar luminosity yields a current mass 
M = 0.89 Mq (Paczyhski 1970), neglecting the possible effects 
of hot-bottom burning. 

Radiative transfer in the inner dust region is properly taken 
into account and the emerging spectrum is used to irradi- 
ate the outer shell. The outer shell is assumed to be opti- 
cally thin. Feed-back from the outer to the inner shell has 
not been taken into account, which is a valid assumption as 
the O-rich outer shell is optically thin at all wavelengths. The 
outer shell contains a mixture of amorphous and crystalline 
silicates. Wo use olivine (Dorschncr, private communication), 
ortho-pyroxene and clino-pyroxene (Koike, private communi- 
cation), and amorphous silicates (Draine & Lee 1984) with 12, 
1, 2 and 85 per cent abundance respectively (by mass), and a 
grain size radius of 0.2 fj,m (0.25 /xm for the amorphous sili- 
cates). The outer radius of the O-rich shell is 4.4 10^'^ cm, but 
is poorly constrained by the model. 

The resulting model fit is shown in Figure 3. Note that our 
solid state model does not yet include all structures seen in the 
spectrum; this will improve with better laboratory data. We 
find mass-loss rates for the C-rich and O-rich shells of 2 10^ ' 
and 8 10~^ Mq/jt respectively, assuming a gas over dust ratio 
of 100 in both shells. The errors in the mass loss are at least 
a factor of two. Reasons for this large uncertainty are, among 



others: (i) the poorly known gas/dust ratio, (ii) the sensitivity 
of M to the grain size distribution function (wc adopt grain 
sizes from 0.008 to 0.1 micron for the amorphous carbon in 
the inner shell and constant values in the outer shell), (Hi) the 
sensitivity of the emission of the O-rich material to the internal 
nebular extinction from the C-rich dust (and the ionized gas). 

The uncertainty in mass loss also refiocts on the accuracy 
of the derived masses, which are ~ 0.01 and ~ 0.37 Mq for 
the inner and outer shell respectively. The O-rich shell mass 
is even more uncertain, as it depends also on the poorly con- 
strained outer boundary radius. The inner radius of the O-rich 
shell corresponds to a dynamical ago of 1050 yr, assuming an 
expansion velocity of 22 km/soc (Acker et al. 1992). Taking 
the O-rich mass at face value, this would imply that this shell 
has been ejected over a time period of approximately 5000 yr. 
This is considerably shorter than the typical inter-puls period. 
The mass loss in the O-rich shell is characteristic for OH/IR 
stars. It may be usefull to search for OH maser emission, as 
one expects that OH emission remains significant until ~ 1000 
to 1500 yr after the end of the AGB. 
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Fig. 3. Model fit to the ISO-SWS spectrum of BD+30 3639. The 
long dashed line is the C-rich dust shell, the short-dashed line rep- 
resents the O-rich detached shell; The dash-dotted line is the sum 
of both components 



4-2. The evolution of [WC] stars 

The SWS spectra shown here demonstrate that both 
BD-f 30 3639 and He 2-113 were O-rich very recently. The age 
of the O-rich AGB remnant (~ 1050 yrs in BD-l-30 3639, and 
possibly less for He 2-113) and the mass of the C-rich enve- 
lope in BD-l-30 3639, suggest that the change to C-rich chem- 
istry was triggered by a PAGB pulse and not by a very late 
thermal pulse. The mass contained in the C-rich dust shell of 
BD+30 3639 is of the same order of magnitude as the sum of 
the H and He layers in a young post- AGB star (Blocker et al. 
1997). This suggests that the PAGB pulse triggered the re- 
moval of the entire outer layer of the star, exposing the He- 
and C-rich inner layers. Evolutionary calculations of Blocker 
& Schonberner (1997) show that in a PAGB thermal pulse, 
no efficient mixing of the H-rich envelope into the He-burning 
layers is expected, and hence no [WC] star would form. Only 
in the case of very late thermal pulses, when the central star is 
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already on the cooling track, mixing of H into the He burning 
layers efficiently removes H, duo to the strongly reduced en- 
velope mass. Since the observations are not consistent with a 
very late thermal pulse for these two nebulae, we suggest that a 
PAGB thermal pulse is capable of mixing H and He-rich layers, 
i.e. contrary to the findings of Blocker & Schonberner (1997). 

Recent calculations of Herwig et al. (1997) suggest that 
a better treatment of convective mixing may result in inter- 
shell abundance patterns consistent with the surface C, O, and 
He abundances of [WC] central stars (see also Leuenhagen & 
Hamann 1997). The latter authors show that late type [WC] 
stars have small amounts of H and N, which may result from 
substantial mixing of H and He-rich layers. It is likely that the 
thermal pulse which caused the change in chemistry is respon- 
sible for this mixing. 

The production of [WC] central stars with 0-rich outer 
shells obviously requires fine-tuning of the timing of the last 
thermal pulse, either at the very end of the AGE or shortly 
thereafter. Zijlstra et al. (1991) discovered OH maser emis- 
sion from IRAS07027-7934, demonstrating that the cool dust 
is 0-rich. Barlow (1998) shows SWS and LWS spectra of 
CPD-56 8032, also with 0-rich dust. This implies that O- 
rich cool dust is not unusual in nebulae with [WC] central 
stars. We note that there is evidence that very late thermal 
pulses can also produce nebulae with [WC] central stars. Pol- 
lacco & Hill (1994) claimed that the [WCll] star 17514-1555 
(PN012.2+04.9) is surrounded by both a very compact nebula 
and an extended, low density nebula, consistent with a very 
late thermal pulse. This is one of the few late-type [WC] stars 
which is relatively faint in the IRAS bands. It is unclear at 
present which mechanism (PAGB thermal pulse with change 
in chemistry or very late thermal pulse) dominates the formar- 
tion rate of PNe with [WC] central stars. 
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